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The hyperfine field from dynamically polarized nuclei in u-GaAs is very spatially inhomogeneous, as the nu¬ 
clear polarization process is most efficient near the randomly-distributed donors. Electrons with polarized spins 
traversing the bulk semiconductor will experience this inhomogeneous hyperfine field as an effective fluctuating 
spin precession rate, and thus the spin polarization of an electron ensemble will relax. A theory of spin relax¬ 
ation based on the theory of random walks is applied to such an ensemble precessing in an oblique magnetic 
field, and the precise form of the (unequal) longitudinal and transverse spin relaxation analytically derived. To 
investigate this mechanism, electrical three-terminal Hanle measurements were performed on epitaxially grown 
Co 2 MnSi/u-GaAs heterostructures fabricated into electrical spin injection devices. The proposed anisotropic 
spin relaxation mechanism is required to satisfactorily describe the Hanle lineshapes when the applied field is 
oriented at large oblique angles. 


Introduction. — The understanding of electrical injection 
and detection of spin in ferromagnetic/semiconductor devices 
has progressed significantly over the past decadeJs^i A key 
obstacle for interpreting spin transport experiments near the 
metal-insulator transition has been the complicating presence 
of dynamically polarized nuclear spins;^“— In the process of 
dynamic nuclear polarization (DNP), the electron spin polar¬ 
ization, maintained out of equilibrium optically or electrically, 
is transferred to the nuclear system over long time scales via 
the hyperfine interaction^—, and can induce nuclear fields up 
to 5.3 T in GaAs. The nature and distribution of the elec¬ 
tronic states controls the properties of the resulting effective 
hyperfine fields from DNP; for instance, electron spins in itin¬ 
erant states interact rapidly with a multitude of nuclei, which 
dilutes the effect and leads to inefficient nuclear polarization. 
Spins situated at impurity sites, however, interact with many 
fewer nuclei, which promotes a more efficient^ DNP. At the 
doping levels examined here, the different donor wave func¬ 
tions overlap often but do not completely fill the bulk crystal, 
which consequently results in a high degree of nuclear field 
inhomogeneity [see Figure [T] ji- i —. Previous descriptions of 
the spin transport dynamics in n-doped semiconductors with 
spin drift-diffusion equations^ii^"— have neglected this essen¬ 
tial inhomogeneity of the nuclear field. 

Here we predict a new anisotropic spin relaxation mech¬ 
anism in semiconductors that occurs when inhomogeneous 
effective magnetic fields are present, such as arise from po¬ 
larized nuclei. Intermediately «-doped GaAs, under the con¬ 
ditions of DNP, offers a testbed for our theory where the 
inhomogeneity manifests itself as a bipartite field with val¬ 
ues Bq or Bq -f Bm with Bq being an applied magnetic field 
and Biy the nuclear field induced by DNP near a donor. We 
demonstrate such a system experimentally and show that mea¬ 
surements of the steady state spin polarization are consistent 
with the devised inhomogeneity-induced anisotropic spin re¬ 
laxation mechanism. 

Theory. — The theoretical description presented herein can 



FIG. 1. (Color online) (a) electron spins (spheres with spin vectors) 
in «-doped GaAs. Bumpy (red) regions depict the presence of donors 
and nuclear fields generated by dynamic nuclear polarization, (b) 
The spin rotation caused by one electron spin entering and departing 
a DNP region, (c) Experimental geometry with 0 being the angle 
between the applied field and sample normal. 


be understood qualitatively by examining Figure [T] wherein 
the inhomogeneity of the magnetic field is represented by ran¬ 
dom placement of distortions (in red) signifying both the pres¬ 
ence of a donor atom and a DNP-induced nuclear field. Elec¬ 
tron spins cross between these donor regions and regions in 
between donors, which lack the DNP-induced nuclear field. 
This transit acts on the spin similarly to an effective fluctuat¬ 
ing Zeeman field. The lower part of Figure [T] shows how the 
inhomogeneity relaxes the spin; in general the nuclear field 
is non-collinear with the applied field, which causes the pre¬ 
cession axis to stochastically modulate when the spin changes 
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field regions. When the transit time between these regions is 
much faster than the change in precession rate experienced 
by the spin upon transit, then the regime of spin relaxation 
corresponds to the motional narrowing regime. We present a 
general form of the theory of spin relaxation in the inhomoge¬ 
neous nuclear field produced by DNP, and specifically within 
the motional narrowing regime we obtain compact analytic re¬ 
sults that can be readily incorporated into spin drift-diffusion 
theories. 

We now present a calculation of the spin relaxation of a 
spin ensemble, 5, ensuing from the aforementioned theory and 
assumptions. In other words we would like to solve for the 
spin relaxation due to the following precession: 


required: 

^=y[Bo+^-bn+^-bnm]xS{t), (3) 

where /(f) is a stochastic function, /(f) is equal to -|-1 (—1) 
for an average time interval \/k„ (1/ko), where 1/ko (l/k„) 
is the average time the spin experiences the field Bq 
before that field changes. We would like to find the dissipative 
effects from the time-dependent field so we will ignore the 
static applied field and average nuclear field: 

^ = lf{t)bnxS{t)=fm-S{t), (4) 


^=y[Bo+Bn{t)]xS{t), (1) 

where the spatial inhomogeneity of the nuclear field is writ¬ 
ten as a time-dependent nuclear field that takes on only two 
possible values of either b„ or 0. Since Bn{t) changes rapidly, 
the first approximation is to replace it with its average value: 
(B„(f)) =Z>„/2 with 


bn — bnuc{I} 


bnucS ■ (Bo + beS) 


{Bo +beS), 


( 2 ) 


where I is the nuclear spin, bnuc is the Overhauser coefficient, 
be is the Knight coefficient, and \/%Bi denotes the strength 
of the random local field. The Knight field allows the nuclear 
field to be non-collinear to the applied held. 

Since the average nuclear held is static, that alone will not 
relax the spin; temporal Huctuations around the average are 

I 


where ii(f) is the skew-symmetric matrix 



/ 0 —cosa sinasinp \ 

— ( cosa 0 —sinacosp 1 , (5) 

y —sinasinp sinacosp 0 / 

where a and p are the spherical coordinates of the nuclear 
held. Depending on the value of /(f), the solution to the pre¬ 
cession equation in between held switchings is 

Sit)=e^-So, Sit)=e-^-So, (6) 

where So is the initial spin vector. 

The time evolution of the spin ensemble can be computed 
by the theory of continuous-time-random-walks . ^ The dif- 
hculty of the theory is reduced since the held modulates be¬ 
tween only two valuesi ^^i^^ The polarization function is a re¬ 
sult of random walks between the two spin environments: 


P(0 = 2 




n{t)+ f -t')e-^^'-''^'Vno{t')e^' dt' + f f O^f - f>^('-'')'To,,(f'- 

Jo Jo Jo 


-I- ... -I- signs of n switched and n o 0 


■So 


r 


where are wait-time distributions to transition between the polarization function in the Laplace domain simplihes to 
state i to state j, and <!>,■ are the survival probabilities in state 


i. Using exponential wait-time distributions leads to: 


P(0 = 2 


Jo 

Jo Jo 

+ ... + signs of n switched and ^ 0 


P{s) = ]- (Ro + R„ + {ko + k„)RoR„) ^ (kok„RoR„)-' • So, 
i >=o 

(8) 

(9) 


with 


R, 


1 


•0(n) 


■So. 


s + ko(„) ± 

which has a Laplace transform equal to 


Utilizing the Laplace transform and its convolution properties. 


f A •-^0 = ^(5)-So. (10) 

s{s + ko+kn)-A. 
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where 


A = -{k„-ko)Cl + SlCl. (11) 


This general expression can be analytically transformed to the 
time domain^ 

We now apply the approximation of fast transitions, 
ko.n >> yBn- To leading order in s, M(s) = [si ~ A/{ko + 
^n)] * which is inverted to be M(f) = gJ^t/ikt+kn) and then 
P{t) = P{t). The next order correction yields^ 


Ht) = ( 


A 

ko “t“ kn 


AA 

{ko + knf 


)pit) 


( 12 ) 


which when written out to second order in becomes 


dP{t) 

dt 


1 

4 ^0 + kn 



y kn kQ ^ 
2 A'o + kn 


xP{t). 


( 13 ) 


Only the first term has the capability to relax the spin ensem¬ 
ble. The second term is a correction to the Larmor precession. 

By combining spin effects such as spin injection, other spin 
relaxation sources, and adding back in the applied and average 
nuclear field in Eq. Q, we can write the following equation 
to encompass the (non-diffusive) spin evolution: 


dP{t) 

dt 


where 


KQ ~\~ kfx 

—P{t) - fibn X (P(f) Xbn) + G, 


X 


1 1 
4 ko + kn 



(14) 


( 15 ) 


Gjjx is the spin generation vector, and is other spin relax¬ 
ation mechanisms which we assume to be isotropic. We have 
simulated the spin evolution with a Monte Carlo approach and 
found agreement with solutions to the differential equation 

Experiment. — To test the theory, we probed the spin polar¬ 
ization P in n-GaAs using electrical Hanle measurements in a 
standard three-terminal (3T) configuration. The sample used 


was an epitaxially grown Co 2 MnSi/«-GaAs (100) heterostruc¬ 
ture. A 2.5 pm thick Si-doped « = 4 x lO'® cm^^ u-GaAs 
channel was grown on an insulating GaAs (100) substrate. To 
thin the naturally occurring Schottky barrier and create a tun¬ 
nel barrier for efficient spin injection^E, a 15 nm n —tran¬ 
sition layer (n = 5 x 10*^ cm^^) was grown followed by a 
18 nm n+ layer. 5 nm of ferromagnetic (FM) Heusler alloy 
Co 2 MnSi was then grown, followed by A1 and Au capping 
layers. 

The structures were patterned into lateral spin injection 
devices^ using standard photolithographic techniques. The 
injection contact was 5 pm x 50 pm. Spin was electrically in¬ 
jected into the n-GaAs channel by imposing a DC current bias 
(800 A/cm^ at 0.51 V) across the FM/n-GaAs interface. For 
the measurements discussed here, the interface was forward- 
biased, so that electrons flowed from the semiconductor into 
the ferromagnet. In a 3T measurement the FM/«-GaAs inter¬ 
face voltage is measured by measuring the voltage with re¬ 
spect to a remote contact outside of the charge current path. 
The spin polarization in the channel directly below the in¬ 
jection contact was probed by measuring the change in the 
3T voltage AV^j upon application of an external out-of-plane 
magnetic field B. This transverse magnetic field serves to pre- 
cess the spins and destroy the spin polarization in the channel 
via the Hanle effecti. 

The influence of DNP on the spin polarization in our sam¬ 
ples is most clearly seen by measuring the 3T Hanle effect 
when the applied field is tilted at small oblique angles 0 away 
from the vertical direction and toward the easy axis of the fer¬ 
romagnetic contact, as shown in Fig. [Uc). The oblique ge¬ 
ometry allows for a significant hyperpolarization of the nuclei 
(Overhauser effect). Satellite peaks are then observed that cor¬ 
respond to fields at which the applied dephasing transverse 
field is partially cancelled by the Overhauser field^. (A less 
prominent satellite peak at very low fields is due to the Knight 
field of the polarized electrons.) The effectiveness in repro¬ 
ducing the oblique 3T Hanle lineshapes therefore serves as a 
test of the validity of the model used to account for the affects 
of DNP. 

Discussion. — Thus far we have only examined the non- 
diffusive dynamics. However the importance of spin diffusion 
on Hanle curves has been well-documented3 In light of the 
theory hitherto presented, we write the following spin diffu¬ 
sion equation: 


^ = y[Bo + X P(f) - -P(t) - y'^bn X iP{t) xbn) + G + Ds^^P + - • VP, (16) 

dt ko -b kn Xi ne 


which is identical to Eq. (fl4l i except for the addition of the 
last two terms which describe spin diffusion and spin drift. 

The physical device geometry was cast into a ID finite- 
element model, where spin may drift and diffuse laterally in 
the sample plane. The simplification to ID is appropriate at 
cryogenic temperatures given the device aspect ratio, where 


the spin diffusion length in GaAs is larger than the chan¬ 
nel thickness. Eq. [16] is iterated forward until steady state 
(^ = 0) is reached. The standard form for the Overhauser 
field^ is used to calculate b„ at each spatial coordinate. Upon 
solving for the steady-state spatially dependent spin polariza¬ 
tion in the channel at each applied field, the 3T Hanle signal 
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(a) 


(b) 




Field (Oe) 


Field (Oe) 


hauser peaks. Table I contains the parameters used to fit the 
3T signal both without and with the anisotropic hyperfine 
terms. Note that the addition of the anisotropic mechanism 
does not change the other fitting parameters significantly, and 
the isotropic lifetime is essentially unchanged. Measurements 
were also taken as a function of injection bias current at fixed 
angle. The fits to the model in this case are comparable to 
those for the angle dependence at fixed bias. Discrepancy be¬ 
tween model and experiment is attributed to a ± 1 ° uncer¬ 
tainty in the angle of field with respect to sample. Addition¬ 
ally the obtained small values for and ko are only on the 
edge of the strong motional narrowing approximation. 


Parameter 

w/o aniso. term 

w/ aniso. term 

Is 

3.3 ns 

3.4 ns 

i^niic 

-1.50 X 10^ Oe 

-1.67 X 10"^ Oe 

be 

-82 Oe 

-73 Oe 

VIBl 

104 Oe 

95 Oe 

ko 

— 

2.1 ns^* 

kfl 

— 

0.45 ns^* 


TABLE I. Fitting parameters for curves in Figure[2l 


FIG. 2. Shown is the oblique-angle 3T Hanle signal measured at 
60 K. The black lines are the experimental data, with a second or¬ 
der (magnetoresistance) background removed and the different an¬ 
gles artificially offset. Shown in gray is the numerical solution to Eq. 
ll6lfor the device geometry, (a) without the anisotropic hyperfine re¬ 
laxation terms and (b) with the anisotropic hyperfine relaxation terms 
included. One set of fitting parameters was used to simultaneously fit 
all the angles. These fitting parameters for both situations are shown 
in Table I. The anisotropic hyperfine terms improved the fit the most 
for larger oblique angles. 

AVar is extracted by projecting the spin polarization at the in¬ 
jector contact Pinj onto the magnetization of the injector fer- 
romagnet M,i i.e AV 37 - 0 = Pi„j ■ M. A single overall scaling 
factor is applied to compare the model to the data. 

In Figure |2] the results of measuring the oblique angle de¬ 
pendence of the 3T Hanle signal at 60 K are shown, along with 
the corresponding fits to the model described above. For com¬ 
parison, the effects of adding the anisotropic hyperfine relax¬ 
ation terms discussed previously are shown side-by-side with 
the fits without the anisotropic hyperfine relaxation terms. In 
both cases, a single set of parameters are used to fit the data 
at all angles. The results show that adding the anisotropic hy¬ 
perfine terms noticeably improve the fitting of the Overhauser 
peak at large oblique angles, for which diffusion alone sys¬ 
tematically overestimates the magnitude and underestimates 
the width of the satellite peak. Without anisotropic relax¬ 
ation, the height and width of the satellite are determined 
only by the spatial variation of the Overhauser field on the 
scale of the electron spin diffusion length. This mechanism 
alone, however, is not sufficient to explain the broadening and 
suppression at larger oblique angles. Inclusion of the addi¬ 
tional smaller length-scale nuclear field inhomogeneity, via 
the anisotropic term, further reduces and broadens the Over- 


Oblique angles larger than ±20° were experimentally in¬ 
accessible due to the switching of the ferromagnetic contact 
when the in-plane component of the field reached the coer¬ 
cive field. Figure 0 shows the solutions of Eq. (fTfil) for two 
larger angles, 30° and 45°. The trend followed at these higher 
angles is similar to what is viewed at the lower ones - the 
anisotropic terms tend to decrease the magnitude of the Over¬ 
hauser peak (black) when compared to their exclusion (red). 
If ferromagnetic contacts with larger coercivity are available, 
a more rigorous test of the predictions of this theory will be 
possible. 



FIG. 3. Large angle solutions to the spin diffusion model, Eq. (ED, 
with and without the anisotropic spin relaxation terms. 

Now we consider how the anisotropic mechanism may also 
be evident in optical spin injection experimentsIn these 
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experiments, the nuclear field is extracted by taking the dif¬ 
ference of the total precession frequency and the precession 
frequency due solely to the applied field.— As we have dis¬ 
cussed here, due to the inherent inhomogeneity of the nuclear 
field, the inferred nuclear field is actually an average nuclear 
field in the probed macroscopic optical spot size. From Eq. 
dl, the inferred nuclear field is then B„ = kob„/{k() +k„) 
which leads to the anisotropic term being 

We predict this term to be observable in time-resolved Faraday 
or Kerr rotation experiments. 


Conclusions. — The influence of DNP on spin evolution in 
semiconductors has been observed for many years. However 
the inherent inhomogeneity of the large nuclear fields has been 
neglected as a spin relaxation process. We have shown that 
the nuclear field inhomogeneity leads to an anisotropic spin 
relaxation mechanism and we have demonstrated that this new 
mechanism can account for the oblique Hanle measurements 
for electrical spin injection into «-GaAs. 
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